Multiple Quantum Oscillations in the de Haas van Alphen Spectra of the Underdoped 

High Temperature Superconductor YBa 2 Cu 3 6 5 
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By improving the experimental conditions and extensive data accumulation, we have achieved 
very high-precision in the measurements of the de Haas- van Alphen effect in the underdoped high- 
temperature superconductor YBa2Cu30e.5. We find that the main oscillation, so far believed to 
be single-frequency, is composed of three closely spaced frequencies. We attribute this to bilayer 
splitting and warping of a single quasi-2D Fermi surface, indicating that c-axis coherence is restored 
at low temperature in underdoped cuprates. Our results do not support the existence of a larger 
frequency of the order of 1650 T reported recently in the same compound [S.E. Sebastian et at, 
Nature 454, 200 (2008)]. 

PACS numbers: 74.25.Bt, 74.25. Ha, 74.72. Bk 



The first unambiguous observation of quantum oscilla- 
tions (QO) in high temperature superconductors (HTSC) 
[l[ has raised much interest concerning the origin of their 
Fermi surface (FS). On the overdoped side, recent quan- 
tum oscillations measurements in Tl 2 Ba 2 Cu06+,s [2[ (Tl- 
2201) have confirmed that the Fermi surface consists of 
a single large cylinder covering 65% of the first Brillouin 
zone (FBZ) as shown by angular magnetoresistance os- 
cillations (AMRO) or angle-resolved photocmission 
spectroscopy (ARPES) 0], in agreement with Hall ef- 
fect data [5j and band structure calculations based on 
local-density approximation (LDA) [fj], 0|. This result 
demonstrates that QO probe directly the FS originating 
from the Cu0 2 planes and that the in-field QO agree 
with zero-field measurements. 

On the underdoped side, both Shubnikov-de Haas 
(SdH) [E 1, HI and de Haas- van Alphen (dHvA) 
[lot 1 1 11 ] oscillations have been observed in underdoped 
YBa2Cu30 !/ (YBCO) pointing to a FS dramatically dif- 
ferent from that of overdoped Tl-2201. Indeed, provided 
that quantum oscillations are interpreted as the conse- 
quence of quantization of closed orbits in a magnetic field 
12], it can be inferred that the FS consists of small pock- 
ets covering only about 2 % of the FBZ jj i, %. One 
fundamental question for our understanding of high tem- 
perature superconductors is what causes the dramatic 
change in FS topology from overdoped to underdoped 
regimes [l3T - fl5l | ? Based on the observation of a negative 
Hall effect at low temperature in underdoped YBCO [l6j |. 
the presence of an electron pocket in the Fermi surface 
has been suggested. It would naturally arise from a re- 
construction of the large hole FS calculated from band 
structure. Different mechanisms have been proposed for 
this FS reconstruction, involving d-DW order (a scenario 
based on a hidden broken symmetry of d^.^-type) [17| . 



incommensurate antiferromagnetic order [ll|, Il8j | , strip e 
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Each scenario predicts a different FS topology and the 
observation of an additional frequency would provide a 
stringent test. In particular, recent dHvA measurements 
flU (performed in the same compound grown at the Uni- 
versity of British Columbia as in ref. [lOj and the present 
study) seems to indicate the presence of a larger pocket, 
which would impose a strong constraint on the ordering 
wave vector for the reconstruction. 

In order to determine whether there are other closed 
sections of the FS, something of fundamental importance 
for clarifying the FS of HTSC in the pseudogap phase, 
we have performed high-precision measurements of the 
de Haas- van Alphen effect in underdoped YBa 2 Cu30 y 
(y=6.51 and y=6.54). We found no evidence of an oscil- 
latory component of high frequency as reported in ref 11 ] 
but the sensitivity of the measurements allow the uncov- 
ering of some important details for the electronic struc- 
ture of underdoped YBCO. In addition to the main fre- 
quency Fi=540 T, two closely spaced satellites F 2 =450 T 
and ^3=630 T have been resolved. We attribute these 
frequencies to a bilayer effect and warping of the FS, 
which proves that quasiparticle motion along the inter- 
plane direction is coherent at low temperature in the un- 
derdoped regime. 

In order to achieve high-precision measurements of 
the de Haas- van Alphen effect, we have averaged several 
pulses (up to 10) for the same experimental conditions, 
as described in more detail in [lOj. We have measured 
two single crystals: YBa 2 Cu306.5i (T c =57.5 K) and 
YBa 2 Cu306.54 (T c =60 K), flux-grown in a non-reactive 
BaZrOs crucible. The dopant oxygen atoms were or- 
dered into the ortho-II superstructure of alternating full 
and empty CuO chains [22j. Typical dimensions of the 
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FIG. 1: (color on line) Oscillatory part of the torque sig- 
nal for underdoped YBa2Cu3 06.54 (upper blue line) and 
YBa2Cu306.5i (lower red line). Solid lines are best fits of 
Eq. [TJ to the data. Four Fourier components are involved in 
the fit (see text). 



samples are 140x140x40 /xm 3 . A careful study of the 
amplitude of the oscillations between the up (50 ms) and 
down (250 ms) field sweeps shows that the sample tem- 
perature changes very little during the pulse. Given an 
effective mass of m* — 1.8 mo (Toj . the upper bound for 
the temperature difference is about 0.2 K at T=0.7 K, 
which excludes any large heating effects during the pulse. 

Fig. Q] displays the field dependence of the oscilla- 
tory torque in two ortho-II samples, obtained after av- 
eraging up to 10 pulses during the down-field sweep at 
T=0.7±0.2 K. A polynomial monotonic background has 
been subtracted from the raw data. While the amplitude 
of the oscillations should grow exponentially with mag- 
netic field for a single frequency, a modulation of the am- 
plitude of the oscillation (beating effect) can be noticed 
in both samples, which is the signature of the existence 
of more than one frequency in the oscillatory spectrum. 
Solid lines in Fig. Q] are fits of the Lifshitz-Kosevich (LK) 
theory, assuming that four frequencies are involved in the 
data, as discussed in more detail below. 

In the framework of the LK theory, the field-dependent 
oscillatory part of the torque, for a multiband 2D Fermi 
surface, can be written: t osc = B Ai sin[27r(-^ — 7,)] 
where Fj is the oscillation frequency linked to a given or- 
bit and 7i is a phase factor. Neglecting magnetic break- 
down and spin damping contributions, the amplitude of 
a given Fourier component Aj oc R^R^i where Rti = 
aTm*/Bsinh[cvTm*/B] and Roi — exp[-aT£)im*/B] arc 
the thermal and Dingle damping factors, respectively. 
Here, a — 27r 2 fcsmo/eS=14.69 Tesla/Kelvin and Tm is 
the Dingle temperature [l2j . In order to restrict the num- 
ber of free parameters in the data analysis, we have con- 
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FIG. 2: (color on line) (a) Oscillatory part of the torque signal 
for YBCOe.54 (symbols). Black solid line is the best fit of 
Eq.[TJto the data with only the main frequency _F\=541 T. (b) 
Residual (fit-data) of the previous fitting (symbols). The line 
shows the best fit of Eq. [TJto the data with two frequencies F\ 
and i<2=453 T. The two other panels show the same procedure 
by taking into account (c) a third frequency F3=637 T and 
(d) a fourth frequency F4 = 1119 T. (e) Residual of the final 
fit. 



sidered the asymptotic form: 

Tosc = ^C l exp(-^)sin[27r(^- - 7*)] (1) 



where B; = a(T+Tm)m* . Strictly speaking, the latter 
assumption is only valid in the high T/B range and only 
affects the amplitude of the signal. The analysis is based 
on discrete Fourier transforms (DFT) using Blackman 
windows, which are known to avoid secondary lobes. 

Fig. displays the field dependence of the oscilla- 
tory torque of YBa 2 Cu3 06.54 (symbols). The solid line is 
a fit to the data assuming only one frequency Fi=541 T 
(solid line), which cannot match the dHvA data over a 
wide field range. This shows that the data cannot be 
accounted for by only one Fourier component. A large 
residual (i.e. the difference between the fit and data) 
can be observed in Fig. f^b) (symbols). The oscillatory 
spectrum still looks complicated, which means that more 
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than one frequency is involved. DFT on this residual re- 
veals the presence of two frequencies close to F\. The 
black line in Fig. ^b) is a fit to the data with a second 
Fourier component ^2=453 T, whose residual is now dis- 
played in Fig. [5Jc) (symbols). The solid line in Fig.[5Jc) 
is a fit with a third frequency ^3=637 T. Finally, the 
residual in Fig. HJd) still contains a small component at 
i<4=1119 T which we assume is the second harmonic of 
F\. The final residual is displayed in Fig.JJJe). The same 
kind of analysis performed on YBa2Cu30g.5i gives sim- 
ilar results. The averaged frequencies deduced from the 
fits of Eq. [1] to the oscillatory torque of the two samples 
are Fi=540±f5 T, F 2 =450±15 T, F 3 =630±40 T and 
F 4 =1I30±20 T. 

The corresponding Fourier analysis of the data for 
YBa2Cu30e.54 and YBa2Cu306.5i are displayed in sym- 
bols in Fig.GUa) and Fig.^b), respectively. A broad peak 
with a maximum around 535 T is observed, in agreement 
with 0, H3]. The DFT of the fit involving only the F x 
component (yellow lines in Fig. [3]) cannot reproduce the 
tail and the small maxima on either side of the broad 
peak. These features are further evidence for the two 
additional frequencies F 2 (green lines in Fig. [3]) and F 3 
(purple lines in Fig. [3]), respectively. The DFT of the fit 
including all the Fourier components successfully repro- 
duces these features, as shown in the inset of Fig. [3ja) 
which displays a zoom of the DFT of the data (symbols) 
and of the fit (black line) for YBa2Cu3C>6.54. The high- 
est frequency detected in both measurements is F4. The 
field dependence of the main peak of the DFT and of the 
oscillation amplitude ('Dingle plot') are displayed in the 
inset of Fig. G^b). DFT of the data (symbols) and DFT 
of the fit (solid lines) are in very good agreement. The 
black lines show the behavior expected if only one fre- 
quency was involved. The strongly non-monotonic field- 
dependence of the frequency of the main peak and the 
apparent non-LK behavior of its amplitude are a signa- 
ture of several Fourier components with frequencies in a 
narrow range. 

Let us now discuss scenarios which may provide an 
explanation for the beating effect observed in QO mea- 
surements. While it cannot be ruled out that additional 
FS pockets account for the observed oscillatory spectra, 
the most natural interpretation is to invoke both a slight 
modulation of the FS sheet along the c-axis (warping) 
and a bilayer splitting effect, the latter being intrinsic 
to YBa2Cu30 y which contains two closely spaced CuC>2 
planes (separated by 3.2 A). The bilayer splitting arises 
from bonding-antibonding coupling of interplane wave 
functions and gives rise to two CuC>2 plane-derived FS 
pockets Q. In band structure calculations for ortho II- 
YBa2Cu 3 06.5 in ref. (23|, the warping and bilayer effects 
led naturally to four frequencies: two come from the bi- 
layer effect and each of them is further split due to warp- 
ing of the FS (providing that c-axis modulation leads to 
two different extremal areas) . Assuming that F4 is a sec- 
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FIG. 3: (color on line) Fourier analysis of the oscillatory 
torque (see Fig.[TJ) for a) YBa2Cu3 06.54 and b) YBa2Cu3C>6.5i 
in the field range 31.4 to 57.4 T. Solid lines show the contribu- 
tion of Fi (yellow), F2 (green) and F3 (purple), respectively. 
The top inset shows a zoom of the DFT of the data (symbols) 
and of the fit including four frequencies (black line). The bot- 
tom inset shows the field dependence of the main peak of the 
DFT and of the oscillation amplitude. Symbols and lines cor- 
respond to DFT of the data and DFT of the fit including four 
frequencies, respectively. 

ond harmonic of Fi (but we can't rule out that F± is a 
second frequency with a small amplitude corresponding 
to the hole pocket in a commensurate FS reconstruction 
[13]), the observation of three frequencies, rather than the 
expected four, can be explained either if two frequencies 
are too close to be resolved or if the amplitude of one 
frequency is too weak. 

There are two hopping terms to consider: t\ comes 
from the hopping through layers in between the CuC>2 
planes and t b ± represents the hopping between the planes 
of the bilayer. In order to estimate t\ which causes 
warping of the FS along the c-axis, we follow the cal- 
culation by Yamaji \l ij and assume the simplest form 
for warping. The quasi-2D dispersion can be written as 
e = h 2 /2m*(kl + ky) — 2t c ± cos(ck z ), where c is the c-axis 
lattice parameter. The c-axis warping leads to two ex- 
tremal FS areas perpendicular to the magnetic field and 
therefore to two close frequencies that yield a beating 
effect. From geometrical considerations, the frequency 
difference is given by AF — Am*t c L /eh. Assuming that 
AF=90 T (F1-F2 or F3-F1), the c-axis hopping term 
tj_ ~f.3 meV, i.e. 15 K. This value is similar to that 
deduced from the beating effect induced by the warp- 
ing of the quasi-2D f3 FS sheet in the quasi-2D oxide 
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Sr 2 Ru0 4 [25 1 . 

In order to deduce a value for the bilayer hopping term 
tj_, it is necessary to rely on band structure calculations. 
A model which assumes a reconstruction of the FS can ac- 
count for the splitting of the frequency linked to the elec- 
tron pocket frequency into two frequencies ^=570 T and 
F=450 T by considering a bilayer coupling t b ± —8 meV 
[26| . For the unreconstructed FS, band structure calcu- 
lations predict that t b ± ~200 meV 0, HI]. This striking 
difference in the value of t b ± is another indication of the 
reconstruction of the large hole FS [39} . 

If warping and bilayer effects are responsible for 
the observation of several frequencies in underdoped 
YBa2Cu3C>6.5, it thus implies that quasiparticle motion 
is coherent along the c-axis at low temperature. While a 
coherent 3D-FS has been established by AMRO experi- 
ments in overdoped Tl-2201 [3J, a crossover in the phase 
diagram between a coherent metal phase (overdoped) and 
an incoherent one (underdoped) has been suggested in 
Bi 2 Sr 2 CaCu 2 8+(5 [27|. Indeed both electrical [H HsJ 
and optical measurements [3CH32| show that in-plane re- 
sponse can be metallic whereas the out-of-plane response 
is insulating-like. The present measurements show that 
at sufficiently low temperature, when superconductivity 
is suppressed by a strong magnetic field, the coherence 
along the c-axis is restored in underdoped and clean ma- 
terial, i.e. it establishes the existence of a coherent three- 
dimensional FS, concomitant with the observation of QO. 

We end by discussing some open questions raised by 
the QO in underdoped YBCO. The first is the apparent 
contradiction with ARPES, which shows the destruction 
of the Fermi surface in the antinodal regions, producing 
a set of disconnected Fermi arcs 



the measurements of the de Haas-van Alphen effect in 
underdoped YBa 2 Cu306.5- We found no evidence of a 
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Assuming that 
ARPES detects only one side of a hole pocket at (tt/2, 

n/2) munis 

is not sufficient to explain any recon- 
struction scenario. Indeed, there is no evidence so far 
in ARPES data of other pockets or regions of the FS 
where coherent quasiparticle exist, except around (tt/2, 
7r/2). In a stripe scenario [19j] ■ the FS contains quasi- 
1D hole-like Fermi surfaces, for which there is no evi- 
dence in ARPES data. One proposed scenario to rec- 
oncile ARPES and QO measurements is a field induced 
state 2(1 21 L a s suggested by recent neutron scattering 
experiment 38j. The second open question concerns the 
larger frequency (Fp ~1650 T), which has been detected 
with a very small amplitude in underdoped YBCO ortho- 
II [HI . The observation of this frequency led the authors 



of ref. IHl to suggest a FS reconstruction of the large hole 
FS into several pockets with different mobilities. In the 
present study, we were not able to detect this high fre- 
quency within a noise level lower than 0.5% of the am- 
plitude of F\ , while the amplitude of this high frequency 
(Fp) is reported to be 3% of the main frequency (Fi) in 
ref. [lH . Further experiments on different samples should 
help clarify this apparent contradiction. 

In summary, we have achieved very high-precision in 



larger pocket as claimed in ref [ll|. A fitting procedure 
based on the Lifshitz-Kosevich theory has revealed that 
the main oscillation is in fact composed of at least three 
closely spaced frequencies, which we attribute to bilayer 
splitting and warping of the FS along the c-axis. This in- 
terpretation implies that quasiparticle transport perpen- 
dicular to the Cu0 2 planes is coherent in the underdoped 
regime at low temperature. 
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